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The Born approximation (Born 1926 Z.Phys.38.802) is a fundamental result in physics, it allows 
the calculation of weak scattering via the Fourier transform of the scattering potential. As was 
done by previous authors (Ge et al 2014 New J. Phys. 16 113048) the Born approximation is 
extended by including in the formula the resonant-states (RSs) of the scatterer. However in this 
study unlike previous studies the included eigen-modes are correctly normalised with dramatic 
positive consequences for the accuracy of the method. The normalisation of the RSs used in the 
previous RSE Born approximation or resonant-state-expansion Born approximation made in Ge et 
al (2014 New J. Phys. 16 113048) has been shown to be numerically unstable in Muljaxov et al 
(2014 arXiv:1409.6877) and by analytics here. The RSs of the system can be calculated using my 
recently discovered RSE perturbation theory for dispersive electrodynamic scatterers (Muljarov et 
al 2010 Europhys. Lett. 92 50010; Boost et al 2012 Phys, Rev. A 89; Boost et al 2014 Phys. Rev. 

A 90 013834) and normalised correctly to appear in the spectral Green’s functions and hence the 
RSE Born approximation via the flux-volume normalisation which I recently rigorously derived in 
Armitage et al (2014 Phys. Rev. A 89), Boost et al (2014 Phys. Rev. A 90 013834)(2016 Phys. Rev. 

A 93 023835). In the case of effectively one-dimensional systems I find an RSE Born approximation 
alternative to the scattering matrix method. 

PACS numbers: 03.50.De, 42.25.-p, 03.65.Nk 

I. INTRODUCTION 

Fundamental to scattering theory, the Born approxi¬ 
mation consists of taking the incident field in place of 
the total field as the driving field at each point inside the 
scattering potential, it was first discovered by Max Born 
and presented in Ref. [l[ . The Born approximation gave 
an expression for the differential scattering cross section 
in terms of the Fourier transform of the scattering po¬ 
tential. The Born approximation is only valid for weak 
scatterers as we will see in the numerical demonstrations. 

In this paper I provide an extension to the Born ap¬ 
proximation which allows an arbitrary number of reso¬ 
nant states (RSs) to be taken into account. I have named 
this extension to the Born approximation the Resonant- 
state-expansion correction to the Born approximation or 
the RSE Born approximation. An almost identical ap¬ 
proach is already available in the literature 0 however 
its derivation differed by including an unstable normal¬ 
isation formula for the RS eigen-modes of the system 
which was then subsequently used to expand Born’s ap¬ 
proximation incorrectly. The normalisation derived in 
Ref.p!^ and used in the previous RSE Born approxima¬ 
tion made in R ef. has been shown to be numerically 
unstable in Ref.[lQ| and shown to be unstable using an¬ 
alytics in Ap^ndix C. Furthermore the numerical study 
made in Ref.[^ only included a single RS in the expan¬ 
sion of the Born approximation, most likely to avoid di¬ 
vergence caused by their incorrect normalisation of the 
RSs. 

Recently there has been developed iii a rigor¬ 
ous perturbation theory called resonant-state expansion 
(RSE) which was then applied to one-dimensional (ID), 


2D and 3D systems 0-0 [HI El which only calculates 
the modes and makes no use of them. The RSE accu¬ 
rately and efficiently calculates RSs of an arbitrary sys¬ 
tem in terms of an expansion of RSs of a simpler, un¬ 
perturbed one. RSs are normalised correctly to appear 
in spectral Green’s functions (GFs) via the flux volume 
normalisation0 and hence the RSE Born approximation. 
In the limit where an infinite number of these resonances 
are included in the RSE Born approximation we will ob¬ 
serve convergence of the method towards the exact solu¬ 
tions. That the RSE can reproduce both the correctly 
normalised RS fields as well as frequencies was demon¬ 
strated in Ref. 0 with the convergence and extrapolation 
algorithm which I contributed to that paper. 

Interestingly the resonant-state-expansion is a near 
identical translation to Electrodynamics of a much ear¬ 
lier theory from Quantum Mechanics by More, Genuw, 
Bang, Gareev, Gizzatkulov and Goncharov Ref.[8|, [Oj]. 
The only difference between the two approaches is the 
choice of RS normalisation method. I am able to show 
in this manuscript that the general normalisation of 
resonant-states which I derived in Ref.0 is the most 
numerically stable available normalisation method. The 
general normalisation which I derived in Ref.0 is based 
on a prototype normalisation which appeared in Ref.0]. 

The concept of RSs was first conceived and used by 
Gamow in 1928 in order to describe mathematically 
the process of radioactive decay, specifically the escape 
from the nuclear potential of an alpha-particle by tun¬ 
nelling. Mathematically this corresponded to solving 
Schrddinger’s equation for outgoing boundary conditions 
(BCs). These states have complex frequency u with 
negative imaginary part meaning their time dependence 
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exp(—zcjt) decays exponentially, thus giving an explana¬ 
tion for the exponential decay law of nuclear physics. The 
consequence of this exponential decay with time is that 
the further from the decaying system at a given instant 
of time the greater the wave amplitude. An intuitive way 
of understanding this divergence of wave amplitude with 
distance is to notice that waves that are further away 
have left the system at an earlier time when less of the 
particle probability density had leaked out. There al¬ 
ready exists numerical techniques for finding eigenmodes 
such as finite element method (FEM) and finite difference 
in time domain (FDTD) method to calculate resonances 
in open cavities. However determining the effect of per¬ 
turbations which break the symmetry presents a signifi¬ 
cant challenge as these popular computational techniques 
need large computational resources to model high quality 
modes. Also these methods generate spurious solutions 
which would damage the accuracy of the RSE Born ap¬ 
proximation if included in the basis. 

The paper is organized as follows, Sec. [IT] outlines the 
derivation of the RSE Born approximation, Sec. [HT] dis¬ 
cusses normalisations of RSs by other authors, Sec. lIVI 
outlines the application of the RSE Born approximation 
to planar slabs. Sec. [Ogives the numerical validation of 
the new method along with a comparison of the alterna¬ 
tive RSE approaches. 


II. DERIVATION OF THE RSE BORN 
APPROXIMATION 

I will in the following section re-derive the method for 
calculating the full GF of an ^en electrodynamic sys¬ 
tem in the same way as Ref.[^ however unlike previ¬ 
ous authors I use the numerically stable normalisation of 
RSs which I derived in Ref.[3, Q. These methods are 
required to calculate transmission and scattering cross- 
section from the dispersive RSE perturbation theory with 
mathematical rigour and accuracy. 

For an electrodynamic system with local frequency de¬ 
pendent dielectric permittivity tensor £fc{r) and perme¬ 
ability /i = 1, where r is the three-dimensional spatial 
position, Maxwell’s wave equation for the electric E(r)- 
field with a current source J(r) oscillating at frequency 
A:, which can be real or complex, is 


- V X V X E{r) + fc2£fc(r)E(r) = ifc—J(r). (1) 

c 

The time-dependent part of the field is given by 
exp(—zcjt) with the complex eigen-frequency cj = ck, 
where c is the speed of light in vacuum. 

The Green’s function (GF) of an open electromagnetic 
system is a tensor which satisfies Maxwell’s wave 
equation Eq. with a delta function source term, 

- Vx Vx Gfc(r,F) + A:^£fc(r)Gfc(r,r') = i^(r-r'), (2) 

where 1 is the unit tensor. Physically, the GF describes 


the response of the system to a point current with fre¬ 
quency uj, i.e. an oscillating dipole. 

The importance of Gfc comes from the fact we can 
see from Eqs. @ that Eqs. o can be solved for E(r) by 
convolution of G^ with the current source J(rO, 

E(r) = y’Gfc(r,r'>jfc^J(r')dr'. (3) 


Inside the system we can use the RSE to calculate the 
GF. In Appendix A I derive for dispersive systems (for 
which I have recently developed a dispersive RSE pertur¬ 
bation theory [l^l a convenient form of the spectral 
GF, valid inside the scatterer only. 


Gfc(r,F) = 

n 


E^(r)(8)E^(F) 

2k{k — kri) 


( 4 ) 


The Eji are RSs of the open optical system and are de¬ 
fined as the eigen-solutions of Maxwell’s wave equation, 

V X V X E„(r) = fc^efc(r)E„(r), (5) 


satisfying the outgoing wave BCs. I have also taken the 
resonator to be embedded in free space (e = 1) without 
loss of generality. Here, is the wave-vector eigen-value 
of the RS numbered by the index n, and En(r) is its elec¬ 
tric field eigen-function. The RSs which are solutions of 
Eq. which are either stationary or dec^in^^in time. 
Modes appearing in Eq. dH) are normalized (7l.ll2l| accord¬ 
ing to the flux-volume normalisation 





d(k‘^ek{r)) 

d{k^) 


E„(r)dr (6) 

k=kn 


-H lim 

k^kn 


En • VE - E • VEr 


Sv 


k‘^-kl 


■dS. 


where the first integral is taken over an arbitrary sim¬ 
ply connected volume V enclosing the inhomogeneity of 
the system and the center of the spherical coordinates 
used, and the second integral is taken over its surface 
Sv- This normalization is required Q for the validity of 
the spectral representation Eq. dJ)- Numerically Eq. dSI) 
has been validated by its use in RSE perturbation the¬ 
ory (S-FtI. [Til , [l ^ . A discussion of the dispersive RSE for 
nano-particles is given in Appendix B. 

The generalisation of my Eq. dH to fcn = 0 modes is 
attributable solely to E. A. Muliarov in Ref. 01 (I derived 
the normalisation proof for Ref. [7| without A: = 0 modes), 
however that part of the proof of the normalisation can 
only be further generalised to dispersive systems using 
the spectral GF Eq. dH) derived in Appendix A. The re¬ 
quired derivation is identical except that it makes use of 
the rigorously derived spectral GF Eq. dJ) instead of the 
identical GF derived in a less mathematically rigorous 
way for non-dispersive systems. This last step is vital for 
the accuracy of the method. Further I note that just as 
I explained in Ref.p!^ must be a symmetric matrix 
or a scalar in order to calculate the dispersion factor as 
shown. 
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Various schemes exist to evaluate the surface inte¬ 
gral limit in Eq. Q such as analytic methods in Ref. [3, 
0? or numerically extending the surface into a non¬ 
reflecting, absorbing, perfectly matched layer where it 
vanishes. 

The derivation of the RSE Born approximation by Ge 
et al. has been made using the normalization intro¬ 
duced by Leung et al. ]^| the limit of infinite volume V 
is taken: 


1 = lim 

V—>CX1 


dk^Skir) 


dik^ 


k=k„ 


KMdr+—f El{r)dS, 
JSv 

( 7 ) 

It was numerically found [13( | that the surface term was 
leading to a stable value of the integral for the relatively 
small volumes available in 2D finite difference in time do¬ 
main (FDTD) calculations. However, it was discovered 
at the time that this was not the case for low-Q modes. 
It was wrongly shown by Muljarov et al [T3| that Eq. Q 
is actually diverging in the limit V —^ oo, and therefore 
the expansion of the Born approximation in and the 
normalization Eq. o are incorrect. In Appendix C I pro¬ 
vided a mathematically rigorous disproof of some of E. 
A. Muljarov’s points and make some correct points about 
the unsuitability of Eq. Q for the RSE perturbation the¬ 
ory myself. Hence although being a cornerstone of the 
scattering theory of open systems the correct expansion 
of the Born approximation in terms of RSs to the exact 
solution was not previously available. 

Analogously to Ref. [3| the derivation of the RSE Born 
approximation of Ge et al. is made but in this case 
using my correct normalisation formula for modes. 

That the E„(r) and k^ can be calculated accurately by 
the RSE perturbation theory and normalised correctly by 
Eq. makes possible the RSE Born approximation. 

The free space GF is now introduced 


- V X V X G{"(r, r') A:^Gf (r, = iS{r - E), (8) 




which has the solution, 


G{‘(r.r') = 


gifc|r—r' I 

47rlr — r' 


( 9 ) 


The systems associated with G^ and G{'® are related 
by the Dyson Equations perturbing back and forth with 
AsfcW = - i 01, 


Combining Eq. (fTUl) and Eq. (fTT|) it is obtained as in 
Ref.i 


Gar,r") = Gf(r.r") 

-ej G7*(r,r')A£.(r')Gf(r'.r")dr' 

+ J j G{''(r,r')A£it(r')Git(r',r"') 

X Ae,(r"')G{“(r"',r")dr"'dr'. (12) 

In order to improve the numerical performance further I 
make a final few steps as in the original Born approxi¬ 
mation [I|, I define unit vector r such that r = rr and 
kg = kv. Then for r » r', 

k\r — E| = fcr ± ks ■ E -K O > 0 (13) 

Therefore substituting Eq. ([1|} and Eq. Q in to Eq. ([T^ 
and using Eq. (fT3|) because both r,r" are far from the 
scatterer we arrive at the RSE Born approximation 


r" 


Gfc(r,E') = - 


47r|r — r" 

gifc(r-+r") 


+ k 


J Ask{r')dr' 

E 


IfiTT^rr" 

V- A„(k.) ® A„(-k'') 

IbTT^rr 


2{k — k„ 


(14) 


or using Eq. 


instead 


r,ik\r—r'‘ 


Gk(ry) = - 
-k^ 


47r|r — r" 

g2fc(r-+r") 




+ k 


IfiTT^rr" 

^ ^.Kr+r") A„(k.) « A„(-k") 

IbTT^rE' ^ 2.kn(k — k^) 


■ (15) 


The vector An is deflned as a Fourier transform of the 
RSs, 


An (ks ) 


J Aek{r')E4r’)dr'. 


(16) 


Gk{r,r'') = 


r") 

(10) 



G{*(r,r"')A£,(r"')G.(r"', 

r")rfr"', 

G,(r"'.r") 

= Gl‘i 

:r'",r") 

(11) 


-k^ 

J Gfc(r"',r')A£fc(r')G{“(r' 

, r")dr', 


I note that the fast Fourier transform method is avail¬ 
able. Furthermore I note that for the inverse scattering 
problem at resonance the inverse Fourier transformation 
is also available. The first two terms in Eq. dTO cor¬ 
respond to the standard Born approximation, the final 
summation term corresponds to the RSE correction to 
the Born approximation. 

A simple corollary of this theory is as follows, we can 
see from the arguments just stated that from Eq.(lO) if 
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y” is inside the resonator and r » r” then 


gife|r-r"| 

47r|r - y'‘ 


A„(k,)«E„(r") 

47rr 2{k — k^) 


(17) 


or using Eq. ()A1[) instead 


Gfc(r,r") = - 


gifc|r-r"| 

47r|r — r'' 


„ikr 




An(ks)(8)E„(r" 


47rr 2kn{k — k„ 


(18) 


similarly from Eq.(ll) if r is inside the resonator and 
r" » r then 


oifc|r-r" I 

47r|r — r" 

^ikr" 

+ k- 


Anr" 

or using Eq. (|A1D instead 

ifc|r-r"| 

Gfc(r,r") = - 


E 


E^(r)® A^(-k ^0 
2(k - kn) 


47r|r — r"\ 

ikr" 




47rr" 2 kn{k — kr, 


(19) 


( 20 ) 


discussing k^ = 0 modes. I also note that k^ = 0 modes 
have only H field or E field component by Maxwell’s 
equations because they are curl free modes, kn = 0 are 
by definition modes which satisfy the condition of being 
curl free. These two points explain why the addition of 
/cn = 0 modes to Sauvan and co-worker’s normalisation 
takes the form it does. Actually the RSE normalisation 
for Hn modes was first shown to me in an email attach¬ 
ment, by E. A. Muljarov several years ago but without 
any derivation and without kn = 0 modes or differential 
dispersion factor included. 

The rigorous derivation of the relationship between 
normalised E-field and normalised H-field can be found 
in appendix D of my PhD thesis [^ . 

Clearly as the perturbation is increased there is a criti¬ 
cal perturbation strength at which the RSE becomes less 
efficient than FDTD and FEM and beyond this point 
one should use the RSE Born approximation with the 
normalisation of Ref.[l3| and FDTD or FEM. 


B. Radiation mode normalisation 

I have recently written a paper on the RSE Born ap¬ 
proximation for waveguides with dispersion fl^ . I found 
that such modes for cylindrical/efFectively-2D waveg¬ 
uides can be normalised by reducing Maxwell’s equation 
to effectively 2D and replacing the operation V x V x with 
a suitable linear operator L invariant along the length of 
the waveguide. A similar approach is found in Ref.[l^ 
and further comparison of the two methods is required. 


other permutations are possible. 


IV. APPLICATION TO PLANAR SYSTEMS 


III. OTHER NORMALISATIONS 


A. Normalisation by Sauvan and co-workers 


The rigorously derived normalisation of Sauvan and 
co-workers that they gave in Ref. [3 


2 = / E„(r)- 

Jv 


d{k^Sk{r)) 


- / H„{r) 

Jv 


d{k‘^) 
d{k‘^ek(r)) 


En{r)dr ( 21 ) 


k=kn 


d(k^ 


Hn(r)dr, 


k=kn 


In this section we discuss the application of the RSE 
Born approximation to exactly solvable ID scattering 
problems in electrodynamics. This is in order to prove 
the converges of the new method to the exact solutions 
available for ID problems in Sec. The dielectric profile 
is described by a scalar frequency independent dielectric 
profile, i.e. €k{z) = le{z), Askiz) = lAs{z). As unper¬ 
turbed system we use a homogeneous planar slab of half 
width a, so that 


{ £5 for \z\ < a , 
1 elsewhere. 


( 22 ) 


requires that the integral V be continued into a perfectly 
matched layer where it is attenuated to zero, thus elim¬ 
inating the need for surface terms in the normalisation. 
As such it is most suitable for FEM and FDTD calcula¬ 
tions. Further I note that just as I explained in Ref. [10] 
£fc(r) must be a symmetric matrix or a scalar in order to 
calculate the dispersion factor as shown. 

I generate this normalisation by combining the RSE 
normalisation for and E^, I note that this alterna¬ 
tive approach to deriving Sauvan’s normalisation was dis¬ 
cussed with E. A. Muljarov at some point but without 


A. Wave equation and normalisation formula in ID 

In this sub-section I consider how Maxwell’s wave 
equation transforms to ID. I also consider how the nor¬ 
malisation formula transforms to ID. 

Maxwell’s wave equations for a planar dielectric struc¬ 
ture with permeability ^ = 1 surrounded by vacuum is 
reduced for ID to the following equation: 

d‘^En{z,t) = £{z)dtEn{z,t) , 


(23) 
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We take the transverse eigen-modes with index n to have 
zero in-plane wave number. The eigen-modes can be fac¬ 
torised as 


= Er^(z)exp{-ickr,t)y (24) 

with time independent part satisfying the wave equation: 

[dl + e{z)kl]E,,{z) = 0, (25) 

The electric field and its first derivative are continuous 
everywhere. Eigenmodes of Maxwell’s wave equation for 
open systems have outgoing boundary conditions. 

In ID non-dispersive systems the RSs En{z) with fre¬ 
quency kn are orthogonal and normalized correctly in ID 
according to [ 3 | 


e[z)En{z)E^{z) dz 

En{-a)Em(-a) E En{a)Em[a) 

i(kn + km) 


= 5r,m , (26) 


where 2 = ±a are the positions of the boundaries of the 
unperturbed system. 


B. Resonant states of the unperturbed slab 


In this sub-section I give the RSs used to calculate the 
RSE Born approximation in Sec. m 

Solving the wave equation Eq. 1)251} for dielectric con¬ 
stant e(z) given by Eq. ([^ . the electric field of RS n, 
normalized according to Eq. (pHll . takes the form Q 


(_l)-A„e-'^-^ z<-a, 

En{z) = { , | 2 | < n, 


An 6 


.‘ik„ z 


n’^ 5 


where 


_ —ifcnd 


An. = 


y/a(es - 1 ) 


with 


2 > o, 


Bn = 


(27) 

(28) 


kn, — 


2a^ 


(ttu. — 2 Iny), n = 0, ±1, ±2, , (29) 


and 


\/^ + 1 
7^’ 


(30) 


the imaginary part of the wave vectors kn are all the 
same. 


C. The form of the RSE Born approximation in 
the one dimensional case 


It is demonstrated in this section that the ID RSE 
Born approximation can be used in conjuncture with the 
RSE perturbation theory (to generate the normalised 
eigen-modes of planar s;^tems with arbitrary dielectric 
profile and dispersion) 1 i, i [H to offer a possible 
alternative to the scattering matrix method of Ref.[T 3 . 
The same method for planar waveguides can be devel¬ 
oped in an analogous w^ except the eigen-modes should 
be calculated as in Ref. p, [I2| • 

In ID the GF Gk(z, z") is the solution of the equation 


ek(z)k^'^Gk{z, z") = S{z- z'‘ 
which from Eq. we can see is given by 

En{z)En{z ) 

^ ^ 2k{k-kn) ■ 


The free space GF is a solution of 

k‘^^Gk(z, z") = S(z - z "), 

and is given by 


Gk[z,z") 



(31) 


(32) 


(33) 


(34) 


Hence in ID the RSE Born approximation is greatly 
simplified to 


Gfc(2,2") = - 


— z) ^ik{z'' —z) pa 

+ -^- 


2ik 


Aek(z)dz 


_ ^ AJk^nik':) 

4 ^ 2{k-kn) ’ 

where An is defined as the Fourier transform, 

An{k,)= [ Aeu(z)Er,[z')dz', (36) 

J —a 


which in the case of a homogeneous slab treated here is 
calculated to be 


i{ks y/^kn) i{ks y/^kn) 

(37)" 

Interestingly in ID we do not require the far field ap¬ 
proximation to make the simplification of the Green’s 
function required to bring the RSE Born approximation 
to the form of Eq. Hence in ID the RSE Born ap¬ 
proximation is valid everywhere outside of the slab and 
not just in the far field. 


An{ks) = Bn As 
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FIG. 1: (Color online) (a) Light transmission through the 
slab. Eq. ((38j. (b) Absolute error in the transmission cal¬ 

culated using the analytic form of r(fc, a) as comparison for 
the numerical values from the RSE Born approximation with 
N = 51, 101, 201,401 as labelled. For further comparison the 
standard Born approximation is also included. 


V. NUMERICAL VALIDATION 

In this section we calculate the ID GF outside of the 
homogeneous slab given by Eq. where = 2.25. We 
do this using the R.SE Born approximation, analytically 
using boundary conditions, and also using the spectral 
GF for comparison. We find that the RSE Born approx¬ 
imation requires less basis states to reach a required ac¬ 
curacy than the spectral GF and unlike the spectral GF 
is convergent outside the system. 

I calculate three types of GF in this section, an an¬ 
alytic GF, the GF of Eq. and the GF of Eq. (i3l^ . 
From these it is possible to use the formula derived in 
Ref. 0, for normal incident and waveguide systems for 
the transmission T{k,z'), 

T{k, z') = \2kG{z', -a; k)f . (38) 

The analytic GF is found by solving Maxwell’s wave 
equation in ID with a source of plane waves while making 
use of Maxwell’s boundary conditions. 

The procedure used to select the basis of RSs for 
the RSE Born approximation calculation is analogous to 
that described in Ref.0 for the RSE perturbation theory. 
Namely, I choose the basis of RSs such that all RSs with 
I^Til < kmax{N) using a maximum wave vector fcmax(N^) 
chosen to select N RSs. 

From Fig. 0 we can see that unlike the standard Born 
approximation the RSE Born approximation is valid over 



zVa 


FIG. 2: Absolute error in transmission calculated using the 
analytic form of T{k,z'), Eq. (I38D . between ^ = —a and z = 
z' as comparison for the RSE Born approximation and the 
spectral GF at fixed ka = 15. A = 201,401 as labelled. 


an arbitrarily wide range of k depending only on the ba¬ 
sis size N used. Furthermore we see that as the basis 
size increases the RSE Born approximation converges to 
the exact solution. The absolute error in the RSE Born 
approximation is approximately reduced by an order of 
magnitude each time the basis size is doubled. Absolute 
errors of 10“^ — 10“^ are seen in the k range shown for 
basis size N = 401. 

From Fig. [2 we can see that unlike the GF calculated 
with the spectral Eq. the RSE Born approximation 
is stable over an arbitrarily wide range of z\ where z’ is 
the coordinate of the point of transmission to, depending 
only on the basis size N used. The transmission cal¬ 
culated via the spectral GF is diverging with distance 
of the point of transmission from the slab, this suggests 
that outside the system the RSE spectral GF is not con¬ 
verging or is divergent. Furthermore we see that the RSE 
Born approximation requires fewer resonant states than 
the spectral GF method in order to produce a required 
absolute error, at all values of z'. Although these points 
were first noted by Ge et al. they were using the un¬ 
stable normalisation leading to an incorrect GF and so 
the results upon which they based their conclusions are 
not reliable. 


VI. SUMMARY 

In this work we have seen the Born approximation 
mathematically rigorously extended to include terms 
which take into account the resonances of the scatter¬ 
ing potential using the exact same method as [2| except 
with correctly normalised modes. Further I have made 
comparisons in ID between scattering calculated with the 
spectral GF and the scattering calculated using the RSE 
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Born approximation. I have demonstrated that once the 
correct normalisation is used in the RSE Born approx¬ 
imation convergences towards the exact solution is ob¬ 
tained. I have found that the RSE Born approximation 
for finding the full GF outside of the system is superior to 
the other spectral GF method considered in terms of con¬ 
vergence and accuracy when the correct normalisation of 
the RSs is used. 

It is demonstrated in this paper that the ID RSE Born 
approximation can be used in conjuncture with the RSE 
perturbation theory (to generate the normalised eigen- 
modes of planar systems with arbitrary dielectric profile) 
[JUSO to offer a possible alternative to the scatter¬ 
ing matrix method of Ref. fl^ . In fact, given the superior 
efficiency of the RSE perturbation theory in comparison 
with FDTD and FEM for weak perturbations demon¬ 
strated in Ref. 01 it is likely that the RSE coupled with 
the RSE Born approximation will be an incredibly pow¬ 
erful scattering theory for weak scatterers. 

I have now derived an analogous theory for general 
wave equations [l^ . 


Appendix A: Derivation of alternative Green’s 
function and completeness 


In order to simplify the RSE Born approximation and 
develop Eq. ([hll we require an appropriate spectral form 
of the GF which is different from the one already proven 
in the literature. To obtain this correct form I start with 
the GF valid inside the scatterer only, which I derived in 

Ref. @.[13, 


Gk{r,T') = 

n 


En(r) 0 En(r^) 


(Al) 


Combining Eq. (1A3D and Eq. (1A4D leads to the closure 
relation 


E„(r') = i'5(r - r'), (A6) 

n 

which expresses the completeness of the RSs, so that any 
function can be written as a superposition of RSs. If in 
the perturbed system some of the series are not conver¬ 
gent or are instead conditionally convergent then we will 
not arrive at the sum rule and completeness, in which 
case I expect that the RSE Born approximation will still 
give convergence to the exact solution but only if a valid 
spectral Green’s function is used, such as Eq. (jAlfl . 


Appendix B: RSE for dispersive systems 


Due to the problems with non-convergence of Schur 
factorisation for the generalised eigen-value problem of 
perturbing nano-spheres dispersively, the RSE in Ref.[Tl| 
might tends to fail for non-symmetric perturbation when 
more than typically 500 basis states are used. This is 
an estimate based on the un-reported RSE failures for 
half and quarter sphere perturbations using the gener¬ 
alised eigen-value problem form of the RSE, tests which 
I carried out for Ref.0. Therefore it is necessary to add 
linear dispersion through a second stage perturbatiom a 
perturbation to the possibly complex conductivity [ll| . 

To make this perturbation consider the problem of a 
perturbation to the conductivity —zA5’(r)/fc^ 


V X V X E^(r) = kl 


- i 


,.A^(r) 


k,. 




(Bl) 


Substituting Eq. (jAl|) in 

-VxVxGfc(r,rO+^^^fc{r)Gfc(r,r') = i(5(r-r^), {A2) 

gives for A: —> oo, 

,(^) ^ (fc + fc„)E„(r) « E„(r') ^ ^ 

ZrZn, 


ekfj, could have in principle any dispersion for which the 
eigen-modes can be normalised, and which becomes non- 
dispersive in the limit of high frequency in order to make 
the sum rule for the GF. 

In this Appendices Greek index letters denote per¬ 
turbed modes and British (English) lower case index let¬ 
ters denote unperturbed modes. 

Since 


since throughout the derivation in this appendix we are 
considering the limit where A: —> oo at which £fc(r) = 
s{v), i.e. the system is non-dispersive at high frequencies. 

Convoluting Eq. (IA3p with arbitrary finite functions 
and assuming the series are convergent we see that since 
A: —> oo we have the sum rule. 


E 


E^(r) (8) E^(r') 
2kn 


= 0 . 


Combining Eq. (lAlD and Eq. (1A4D yields 


(A4) 


I Gfc^(r,r')A<T(rOE^(r')dr', (B2) 


where (r, r') is given by Eq. (lAfjp and by Eq. 

E^r) = ^6„^E„(r). (B3) 


where in Eq. ()B3[) E„ and k^ correspond to the unper¬ 
turbed modes of Eq. ([fi]), then following the derivation 
method of Ref. 0 0 


Gfc(r.r') = E 

n 


E^(r) 0 E^(r^) 
2A:(A: — k^) 


(AS) 


^^nfikfj, = E (^®TH2 “1“ ^^naka) ^afi ; (®4) 

a 










where 


S„a = j E„(r) . A,T(r)E„(r)dr, (B5) 

which can be solved for the eigen-modes and eigen¬ 
values kfj_ of the perturbed problem. 

For examples of fitting dispersion linear in wavelength 
to the dispersion of real materials for the purposes of RSE 
perturbation theory please see my Ref.(ll| where a lin¬ 
ear dispersive RSE is presented in terms of a generalised 
eigen-value problem. 

The perturbation to can also be non-dispersive 
without resorting to generalised eigen-value problems, as 
treated in my Ref.Q]. To elaborate further on this point 
consider the problem of the non-dispersive perturbation 
Ae(r) 

V X V X (r) = fcjl [Sk^ (r) + A£(r)] E^ (r). (B 6 ) 

Again Sk is dispersive as in Eq. Eq. (IB 6 II is solved 

by 0 



where 

= j E„(r) ■ Ae(r)E„(r)dr (B 8 ) 

and bnfj,y/k^ = Please note that it is very impor¬ 
tant to be consistent with the signs of \/^ in the matrix 
elements of Eq. (1B7I1 . 

Using the linear eigen-value approach outlined here it 
might be possible to treat an unperturbed Drude-Lorentz 
gold sphere with a non-dispersive shell, and perturb away 
the non-dispersive shell leaving in its place biological par¬ 
ticles to be sensed as a perturbation. All perturbations 
must be within the boundaries of the unperturbed system 
due to convergence of the GF, see Fig. El 

For a discussion of the eigen-functions of Maxwell’s 
equations in spherical coordinates please see [ 2 ^. 


Appendix C: Kristensen normalisation 

In order for the normalisation of Kristensen et al to be 
correct it must be consistent with my Eq. (I^, specifically 
the surface term in Eq. Q must be mathematically 

equivalent with Kristensen’s surface term . Hence it 
should be that where from Q we have, 


(Cl) 

= (C2) 

V ->00 Ikn J 


<t> dS 


Js„ 


„ d 9E„ 


dr 


and from Eq. o we have 


However considering the RSE Born approximation in 3D, 
specifically to ensure outgoing BCs Eq. CZD in particular 
driven by (convoluted with) a current vanishing propor¬ 
tional to {k — kn) as fc —>• so E(r) E„(r), we know 
that 


lim E„{r) =£{«,</,)—. (C3) 

r—>^00 T 

where r = (r, 9, (j)) in spherical polar coordinates and 
f(^,(/)) oc Aji(ffc„), then substituting Eq. (1C3II into 
Eq. and Eq. and equating and gives, 

lim = lim ^^ </ f^(0,0)^5', (C4) 

r—>^00 r—^00 zkjiT^ JSv 

a logically valid statement, therefore, when 

r —¥ 00 and so the normalisation of Kristensen et al is not 
wrong as it is stated. In the RSE perturbation theory 
letting r —> 00 in the normalisation of perturbed modes 
introduces huge errors because of the blow up of the RS 
mode fields far from the system causing blow up of er¬ 
ror. By Eq. (|C4|1 . as r grows one is essentially subtracting 
an exponentially growing surface term from an exponen¬ 
tially growing volume term to get the constant 1 - 1 - Sk^,o 
for normalisation, this leads to large numerical errors for 
low-Q (leaky) modes [3. An inherent source of insta¬ 
bility is remaining dependence of f(^, 0 ) on r due to the 
use of finite r. 

In ID Kristensen et al normalisation is actually correct 
for any finite r that includes the system inhomogeneity 

i- 

The remaining problems with the Kristensen et al nor¬ 
malisation is that it is missing k = 0 modes, therefore 
it is incomplete, and hence incorrect. Also it does not 
have the conditions on which should be the same as 
for my normalisation. That the RSs can be written in 
the form of Eq. (|C3D aids the solution of the inverse scat¬ 
tering problem [ 12 j. 

As an aside, since 


f(0,0) = (C5) 

47r 

it is theoretically possible to make calculations of the po¬ 
tential from the emission (decay) via fast inverse Fourier 
transform methods upon the set of An(rfcn) especially 
if the potentials of interest are rotating about a fixed 
axis so we know their orientation to some extent such as 
occurs for decaying magnetic nuclei as part of a non¬ 
magnetic crystalline compound placed inside a NMR 
(nuclear-magnetic-resonance) machine. Because k^ are 
discrete values these inverse Fourier methods might have 
to be used self-consistently in conjuncture with the RSE 
perturbation theory and the values of k^. This is a highly 
speculative aside and might be a possible topic for futme 
research. 

As another aside, close to a sharp resonance k ^ kn 
scattering from the potential is dominated by a single 
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resonance and so the scattered E-field is ap¬ 

proximately 


and the arguments with regards to normalisation in this 
Appendix C still hold except that now 


lim 

r—¥oo 


pikr 

(r, fc) CA„(rfc)-, 

r 


{C6) 


again A„(rA:) can be partially inverse Fourier trans¬ 
formed with respect to angle to find information about 
the internal structure of the potential. C is some con¬ 
stant. This argument assumed elastic scattering, which 
is a valid assumption even for such things as neutron- 
nucleus scattering provided that the neutron energies are 
high enough, at low energies inelastic scattering causes 
deviations from the RSE Born approximation model, 
these deviations in scattering caused by inelasticity give 
the fission or absorption cross-sections. An RSE for 
Schrodinger’s equation is given in [I^ . 

For 2D systems, by following similar arguments as here 
we arrive at (in the notation of 0 ) 

lim E„(r) = , (C7) 

P^oo ^pTT 


1 + — J ■ 


d{uj^e 




E„.dr 


(C8) 


R 


27 r 


de 


Erj • —-h pEr. 


dp 


dp‘^ 




Jp=i2 
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